Genome-wide diversity, population structure and signatures of inbreeding in the African buffalo in Mozambique.

The African buffalo, Syncerus caffer, is a key species in African ecosystems. Like other large herbivores, it plays a fundamental role in its habitat acting as an ecosystem engineer. Over the last few centuries, African buffalo populations have declined because of range contraction and demographic decline caused by direct or indirect human activities. In Mozambique, historically home to large buffalo herds, the combined effect of colonialism and subsequent civil wars has created a critical situation that urgently needs to be addressed. In this study, we focused on the analysis of genetic diversity of Syncerus caffer caffer populations from six areas of Mozambique. Using genome-wide SNPs obtained from ddRAD sequencing, we examined the population structure across the country, estimated gene flow between areas under conservation management, including national reserves, and assessed the inbreeding coefficients. Our results indicate that all studied populations of Syncerus caffer caffer are genetically depauperate, with a high level of inbreeding. Moreover, buffaloes in Mozambique present a significant population differentiation between southern and central areas. We found an unexpected genotype in the Gorongosa National Park, where buffaloes experienced a dramatic population size reduction, that shares a common ancestry with southern populations of Catuane and Namaacha. This could suggest the past occurrence of a connection between southern and central Mozambique and that the observed population structuring could reflect recent events of anthropogenic origin. All the populations analysed showed high levels of homozygosity, likely due to extensive inbreeding over the last few decades, which could have increased the frequency of recessive deleterious alleles. Improving the resilience of Syncerus caffer caffer in Mozambique is essential for preserving the ecosystem integrity. The most viable approach appears to be facilitating translocations and re-establishing connectivity between isolated herds. However, our results also highlight the importance of assessing intraspecific genetic diversity when considering interventions aimed at enhancing population viability such as selecting suitable source populations.

Comparative genomics provides insights into molecular adaptation to hypermetamorphosis and cantharidin metabolism in blister beetles (Coleoptera: Meloidae).

Blister beetles (Coleoptera: Meloidae) are currently subdivided into three subfamilies: Eleticinae (a basal group), Nemognathinae, and Meloinae. These are all characterized by the endogenous production of the defensive terpene cantharidin (CA), whereas the two most derived subfamilies show a hypermetamorphic larval development. Here, we provide novel draft genome assemblies of five species sampled across the three blister beetle subfamilies (Iselma pallidipennis, Stenodera caucasica, Zonitis immaculata, Lydus trimaculatus, and Mylabris variabilis) and performed a comparative analysis with other available Meloidae genomes and the closely-related canthariphilous species (Pyrochroa serraticornis) to disclose adaptations at a molecular level. Our results highlighted the expansion and selection of genes potentially responsible for CA production and metabolism, as well as its mobilization and vesicular compartmentalization. Furthermore, we observed adaptive selection patterns and gain of genes devoted to epigenetic regulation, development, and morphogenesis, possibly related to hypermetamorphosis. We hypothesize that most genetic adaptations occurred to support both CA biosynthesis and hypermetamorphosis, two crucial aspects of Meloidae biology that likely contributed to their evolutionary success.

Limited introgression from non-native commercial strains and signatures of adaptation in the key pollinator Bombus terrestris.

Insect pollination is fundamental for natural ecosystems and agricultural crops. The bumblebee species Bombus terrestris has become a popular choice for commercial crop pollination worldwide due to its effectiveness and ease of mass rearing. Bumblebee colonies are mass produced for the pollination of more than 20 crops and imported into over 50 countries including countries outside their native ranges, and the risk of invasion by commercial non-native bumblebees is considered an emerging issue for global conservation and biological diversity. Here, we use genome-wide data from seven wild populations close to and far from farms using commercial colonies, as well as commercial populations, to investigate the implications of utilizing commercial bumblebee subspecies in the UK. We find evidence for generally low levels of introgression between commercial and wild bees, with higher admixture proportions in the bees occurring close to farms. We identify genomic regions putatively involved in local and global adaptation, and genes in locally adaptive regions were found to be enriched for functions related to taste receptor activity, oxidoreductase activity, fatty acid and lipid biosynthetic processes. Despite more than 30 years of bumblebee colony importation into the UK, we observe low impact on the genetic integrity of local B. terrestris populations, but we highlight that even limited introgression might negatively affect locally adapted populations.

Comparative ontogenetic and transcriptomic analyses shed light on color pattern divergence in cichlid fishes.

Stripe patterns are a striking example for a repeatedly evolved color pattern. In the African adaptive radiations of cichlid fishes, stripes evolved several times independently. Previously, it has been suggested that regulatory evolution of a single gene, agouti-related-peptide 2 (agrp2), explains the evolutionary lability of this trait. Here, using a comparative transcriptomic approach, we performed comparisons between (adult) striped and nonstriped cichlid fishes of representatives of Lake Victoria and the two major clades of Lake Malawi (mbuna and non-mbuna lineage). We identify agrp2 to be differentially expressed across all pairwise comparisons, reaffirming its association with stripe pattern divergence. We therefore also provide evidence that agrp2 is associated with the loss of the nonstereotypic oblique stripe of Mylochromis mola. Complementary ontogenetic data give insights into the development of stripe patterns as well as vertical bar patterns that both develop postembryonically. Lastly, using the Lake Victoria species pair Haplochromis sauvagei and Pundamilia nyererei, we investigated the differences between melanic and non-melanic regions to identify additional genes that contribute to the formation of stripes. Expression differences-that most importantly also do not include agrp2-are surprisingly small. This suggests, at least in this species pair, that the stripe phenotype might be caused by a combination of more subtle transcriptomic differences or cellular changes without transcriptional correlates. In summary, our comprehensive analysis highlights the ontogenetic and adult transcriptomic differences between cichlids with different color patterns and serves as a basis for further investigation of the mechanistic underpinnings of their diversification.

Comparative Genomics of the Waterfowl Innate Immune System.

Animal species differ considerably in their ability to fight off infections. Finding the genetic basis of these differences is not easy, as the immune response is comprised of a complex network of proteins that interact with one another to defend the body against infection. Here, we used population- and comparative genomics to study the evolutionary forces acting on the innate immune system in natural hosts of the avian influenza virus (AIV). For this purpose, we used a combination of hybrid capture, next- generation sequencing and published genomes to examine genetic diversity, divergence, and signatures of selection in 127 innate immune genes at a micro- and macroevolutionary time scale in 26 species of waterfowl. We show across multiple immune pathways (AIV-, toll-like-, and RIG-I -like receptors signalling pathways) that genes involved genes in pathogen detection (i.e., toll-like receptors) and direct pathogen inhibition (i.e., antimicrobial peptides and interferon-stimulated genes), as well as host proteins targeted by viral antagonist proteins (i.e., mitochondrial antiviral-signaling protein, [MAVS]) are more likely to be polymorphic, genetically divergent, and under positive selection than other innate immune genes. Our results demonstrate that selective forces vary across innate immune signaling signalling pathways in waterfowl, and we present candidate genes that may contribute to differences in susceptibility and resistance to infectious diseases in wild birds, and that may be manipulated by viruses. Our findings improve our understanding of the interplay between host genetics and pathogens, and offer the opportunity for new insights into pathogenesis and potential drug targets.

Bacterial microcompartments for isethionate desulfonation in the taurine-degrading human-gut bacterium Bilophila wadsworthia.

BACKGROUND: Bilophila wadsworthia, a strictly anaerobic, sulfite-reducing bacterium and common member of the human gut microbiota, has been associated with diseases such as appendicitis and colitis. It is specialized on organosulfonate respiration for energy conservation, i.e., utilization of dietary and host-derived organosulfonates, such as taurine (2-aminoethansulfonate), as sulfite donors for sulfite respiration, producing hydrogen sulfide (H2S), an important intestinal metabolite that may have beneficial as well as detrimental effects on the colonic environment. Its taurine desulfonation pathway involves the glycyl radical enzyme (GRE) isethionate sulfite-lyase (IslAB), which cleaves isethionate (2-hydroxyethanesulfonate) into acetaldehyde and sulfite. RESULTS: We demonstrate that taurine metabolism in B. wadsworthia 3.1.6 involves bacterial microcompartments (BMCs). First, we confirmed taurine-inducible production of BMCs by proteomic, transcriptomic and ultra-thin sectioning and electron-microscopical analyses. Then, we isolated BMCs from taurine-grown cells by density-gradient ultracentrifugation and analyzed their composition by proteomics as well as by enzyme assays, which suggested that the GRE IslAB and acetaldehyde dehydrogenase are located inside of the BMCs. Finally, we are discussing the recycling of cofactors in the IslAB-BMCs and a potential shuttling of electrons across the BMC shell by a potential iron-sulfur (FeS) cluster-containing shell protein identified by sequence analysis. CONCLUSIONS: We characterized a novel subclass of BMCs and broadened the spectrum of reactions known to take place enclosed in BMCs, which is of biotechnological interest. We also provided more details on the energy metabolism of the opportunistic pathobiont B. wadsworthia and on microbial H2S production in the human gut.

Phosphitispora fastidiosa gen. nov. sp. nov., a new dissimilatory phosphite-oxidizing anaerobic bacterium isolated from anaerobic sewage sludge.

A new strictly anaerobic bacterium, strain DYL19T, was enriched and isolated with phosphite as the sole electron donor and CO2 as a single carbon source and electron acceptor from anaerobic sewage sludge sampled at a sewage treatment plant in Constance, Germany. It is a Gram-positive, spore-forming, slightly curved, rod-shaped bacterium which oxidizes phosphite to phosphate while reducing CO2 to biomass and small amounts of acetate. Optimal growth is observed at 30 °C, pH 7.2, with a doubling time of 3 days. Beyond phosphite, no further inorganic or organic electron donor can be used, and no other electron acceptor than CO2 is reduced. Sulphate inhibits growth with phosphite and CO2. The G+C content is 45.95 mol%, and dimethylmenaquinone-7 is the only quinone detectable in the cells. On the basis of 16S rRNA gene sequence analysis and other chemotaxonomic properties, strain DYL19T is described as the type strain of a new genus and species, Phosphitispora fastidiosa gen. nov., sp. nov.

Rapid adaptive radiation in a hillstream cyprinid fish in the East African White Nile River basin.

Adaptive radiation of freshwater fishes was long thought to be possible only in lacustrine environments. Recently, several studies have shown that riverine and stream environments also provide the ecological opportunity for adaptive radiation. In this study, we report on a riverine adaptive radiation of six ecomorphs of cyprinid hillstream fishes of the genus Garra in a river located in the Ethiopian Highlands in East Africa. Garra are predominantly highly specialized algae-scrapers with a wide distribution ranging from Southeast Asia to West Africa. However, adaptive phenotypic diversification in mouth type, sucking disc morphology, gut length and body shape have probably been found among these ecomorphs in a single Ethiopian river. Moreover, we found two novel phenotypes of Garra ("thick-lipped" and "predatory") that had not been discovered before in this species-rich genus (>160 species). Mitochondrial and genome-wide data suggest monophyletic, intrabasin evolution of Garra phenotypic diversity with signatures of gene flow from other local populations. Although sympatric ecomorphs are genetically distinct and can be considered to being young species as suggested by genome-wide single nucleotide polymorphism data, mitochondrial DNA was unable to identify any genetic structure suggesting recent and rapid speciation events. Some data suggest a hybrid origin of the novel "thick-lipped" ecomorph. Here we highlight how, driven by ecological opportunity, an ancestral trophically highly specialized lineage is likely to have rapidly radiated in a riverine environment promoted by the evolution of novel feeding strategies.

Pelorhabdus rhamnosifermentans gen. nov., sp. nov., a strictly anaerobic rhamnose degrader from freshwater lake sediment.

A rhamnose-degrading bacterium, strain BoRhaAT, was isolated from profundal sediment of Lake Constance in agar dilution series with l-rhamnose as substrate and with a background lawn of Methanospirillum hungatei. The isolated strain was a motile rod that stained Gram positive. Growth was observed within a pH range of 4.0-7.5 and a temperature range of 15-30°C. Fermentation products of rhamnose or glucose were acetate, propionate, ethanol, butyrate, and 1-propanol. The G+C content was 40.6% G+C. The dominant fatty acids are C16:1ω9c, i-C13:03OH, C16:0 and C17:1ω8c with 8-21% relative abundance. Polar lipids were glycolipids, phosphatidylethanolamine, phosphoaminolipid and other lipids, of which phosphatidylethanolamine was most abundant. The sequence of the 16S rRNA gene of the new isolate matches the sequence of its closest relative Anaerosporomusa subterranea to 92.4%. A comparison of the genome with this strain showed 60.2% genome-wide average amino acid identity (AAI), comparisons with other type strains showed a maximum of 62.7% AAI. Thus, the definition of a new genus is justified for which we propose the name Pelorhabdus. For strain BoRhaAT, we propose the name Pelorhabdus rhamnosifermentans gen. nov., sp. nov., with strain BoRhaAT (DSM 111565T = JCM 39158T) as the type strain.

Capturing the rapidly evolving study of adaptation.

Research on the genomics of adaptation is rapidly changing. In the last few decades, progress in this area has been driven by methodological advances, not only in the way increasingly large amounts of molecular data are generated (e.g. with high-throughput sequencing), but also in the way these data are analysed. This includes a growing appreciation and quantitative treatment of covariation among units within the same data type (e.g. genes) or across data types (e.g. genes and phenotypes). The development and adoption of more and more integrative tools have resulted in richer and more interesting empirical work. This special issue - comprising methodological, empirical, and review papers - aims to capture a 'snapshot' of this rapidly evolving field. We discuss in particular three important themes in the study of adaptation: the genetic architecture of adaptive variation, protein-coding and regulatory changes, and parallel evolution. We highlight how more traditional key themes in the study of genetic architecture (e.g. the number of loci underlying adaptive traits and the distribution of their effects) are now being complemented by other factors (e.g. how patterns of linkage and number of loci interact to affect the ability to adapt). Similarly, apart from addressing the relative importance of protein-coding and regulatory changes, we now have the tools to look in-depth at specific types of regulatory variation to gain a clearer picture of regulatory networks. Finally, parallel evolution has always been central to the study of adaptation, but now we are often able to address the question of whether - and to what extent - parallelism at the organismal or phenotypic level is matched by parallelism at the genetic level. Perhaps most importantly, we can now determine what mechanisms are driving parallelism (or lack thereof) across levels of biological organization. All these recent methodological developments open up new directions for future studies of adaptive changes across traits, levels of biological organization, demographic contexts and time scales.

Towards complete and error-free genome assemblies of all vertebrate species.

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

Giant lungfish genome elucidates the conquest of land by vertebrates.

Lungfishes belong to lobe-fined fish (Sarcopterygii) that, in the Devonian period, 'conquered' the land and ultimately gave rise to all land vertebrates, including humans1-3. Here we determine the chromosome-quality genome of the Australian lungfish (Neoceratodus forsteri), which is known to have the largest genome of any animal. The vast size of this genome, which is about 14× larger than that of humans, is attributable mostly to huge intergenic regions and introns with high repeat content (around 90%), the components of which resemble those of tetrapods (comprising mainly long interspersed nuclear elements) more than they do those of ray-finned fish. The lungfish genome continues to expand independently (its transposable elements are still active), through mechanisms different to those of the enormous genomes of salamanders. The 17 fully assembled lungfish macrochromosomes maintain synteny to other vertebrate chromosomes, and all microchromosomes maintain conserved ancient homology with the ancestral vertebrate karyotype. Our phylogenomic analyses confirm previous reports that lungfish occupy a key evolutionary position as the closest living relatives to tetrapods4,5, underscoring the importance of lungfish for understanding innovations associated with terrestrialization. Lungfish preadaptations to living on land include the gain of limb-like expression in developmental genes such as hoxc13 and sall1 in their lobed fins. Increased rates of evolution and the duplication of genes associated with obligate air-breathing, such as lung surfactants and the expansion of odorant receptor gene families (which encode proteins involved in detecting airborne odours), contribute to the tetrapod-like biology of lungfishes. These findings advance our understanding of this major transition during vertebrate evolution.

The comparative genomic landscape of adaptive radiation in crater lake cichlid fishes.

Factors ranging from ecological opportunity to genome composition might explain why only some lineages form adaptive radiations. While being rare, particular systems can provide natural experiments within an identical ecological setting where species numbers and phenotypic divergence in two closely related lineages are notably different. We investigated one such natural experiment using two de novo assembled and 40 resequenced genomes and asked why two closely related Neotropical cichlid fish lineages, the Amphilophus citrinellus species complex (Midas cichlids; radiating) and Archocentrus centrarchus (Flyer cichlid; nonradiating), have resulted in such disparate evolutionary outcomes. Although both lineages inhabit many of the same Nicaraguan lakes, whole-genome inferred demography suggests that priority effects are not likely to be the cause of the dissimilarities. Also, genome-wide levels of selection, transposable element dynamics, gene family expansion, major chromosomal rearrangements and the number of genes under positive selection were not markedly different between the two lineages. To more finely investigate particular subsets of the genome that have undergone adaptive divergence in Midas cichlids, we also examined if there was evidence for 'molecular pre-adaptation' in regions identified by QTL mapping of repeatedly diverging adaptive traits. Although most of our analyses failed to pinpoint substantial genomic differences, we did identify functional categories containing many genes under positive selection that provide candidates for future studies on the propensity of Midas cichlids to radiate. Our results point to a disproportionate role of local, rather than genome-wide factors underlying the propensity for these cichlid fishes to adaptively radiate.

Does divergence from normal patterns of integration increase as chromosomal fusions increase in number? A test on a house mouse hybrid zone.

Chromosomal evolution is widely considered an important driver of speciation because it can promote the establishment of reproductive barriers. Karyotypic reorganization is also expected to affect the mean phenotype, as well as its development and patterns of phenotypic integration, through processes such as variation in genetic linkage between quantitative trait loci or between regulatory regions and their targets. Here we explore the relationship between chromosomal evolution and phenotypic integration by analyzing a well-known house mouse parapatric contact zone between a highly derived Robertsonian (Rb) race (2n = 22) and populations with standard karyotype (2n = 40). Populations with hybrid karyotypes are scattered throughout the hybrid zone connecting the two parental races. Using mandible shape data and geometric morphometrics, we test the hypothesis that patterns of integration progressively diverge from the "normal" integration pattern observed in the standard race as they accumulate Rb fusions. We find that the main pattern of integration observed between the posterior and anterior part of the mandible can be largely attributed to allometry. We find no support for a gradual increase in divergence from normal patterns of integration as fusions accumulate. Surprisingly, however, we find that the derived Rb race (2n = 22) has a distinct allometric trajectory compared with the standard race. Our results suggest that either individual fusions disproportionately affect patterns of integration or that there are mechanisms which "purge" extreme variants in hybrids (e.g. reduced fitness of hybrid shape).

Contrasting signatures of genomic divergence during sympatric speciation.

The transition from 'well-marked varieties' of a single species into 'well-defined species'-especially in the absence of geographic barriers to gene flow (sympatric speciation)-has puzzled evolutionary biologists ever since Darwin1,2. Gene flow counteracts the buildup of genome-wide differentiation, which is a hallmark of speciation and increases the likelihood of the evolution of irreversible reproductive barriers (incompatibilities) that complete the speciation process3. Theory predicts that the genetic architecture of divergently selected traits can influence whether sympatric speciation occurs4, but empirical tests of this theory are scant because comprehensive data are difficult to collect and synthesize across species, owing to their unique biologies and evolutionary histories5. Here, within a young species complex of neotropical cichlid fishes (Amphilophus spp.), we analysed genomic divergence among populations and species. By generating a new genome assembly and re-sequencing 453 genomes, we uncovered the genetic architecture of traits that have been suggested to be important for divergence. Species that differ in monogenic or oligogenic traits that affect ecological performance and/or mate choice show remarkably localized genomic differentiation. By contrast, differentiation among species that have diverged in polygenic traits is genomically widespread and much higher overall, consistent with the evolution of effective and stable genome-wide barriers to gene flow. Thus, we conclude that simple trait architectures are not always as conducive to speciation with gene flow as previously suggested, whereas polygenic architectures can promote rapid and stable speciation in sympatry.

The mole genome reveals regulatory rearrangements associated with adaptive intersexuality.

Linking genomic variation to phenotypical traits remains a major challenge in evolutionary genetics. In this study, we use phylogenomic strategies to investigate a distinctive trait among mammals: the development of masculinizing ovotestes in female moles. By combining a chromosome-scale genome assembly of the Iberian mole, Talpa occidentalis, with transcriptomic, epigenetic, and chromatin interaction datasets, we identify rearrangements altering the regulatory landscape of genes with distinct gonadal expression patterns. These include a tandem triplication involving CYP17A1, a gene controlling androgen synthesis, and an intrachromosomal inversion involving the pro-testicular growth factor gene FGF9, which is heterochronically expressed in mole ovotestes. Transgenic mice with a knock-in mole CYP17A1 enhancer or overexpressing FGF9 showed phenotypes recapitulating mole sexual features. Our results highlight how integrative genomic approaches can reveal the phenotypic impact of noncoding sequence changes.

Environmental and Intestinal Phylum Firmicutes Bacteria Metabolize the Plant Sugar Sulfoquinovose via a 6-Deoxy-6-sulfofructose Transaldolase Pathway.

Bacterial degradation of the sugar sulfoquinovose (SQ, 6-deoxy-6-sulfoglucose) produced by plants, algae, and cyanobacteria, is an important component of the biogeochemical carbon and sulfur cycles. Here, we reveal a third biochemical pathway for primary SQ degradation in an aerobic Bacillus aryabhattai strain. An isomerase converts SQ to 6-deoxy-6-sulfofructose (SF). A novel transaldolase enzyme cleaves the SF to 3-sulfolactaldehyde (SLA), while the non-sulfonated C3-(glycerone)-moiety is transferred to an acceptor molecule, glyceraldehyde phosphate (GAP), yielding fructose-6-phosphate (F6P). Intestinal anaerobic bacteria such as Enterococcus gilvus, Clostridium symbiosum, and Eubacterium rectale strains also express transaldolase pathway gene clusters during fermentative growth with SQ. The now three known biochemical strategies for SQ catabolism reflect adaptations to the aerobic or anaerobic lifestyle of the different bacteria. The occurrence of these pathways in intestinal (family) Enterobacteriaceae and (phylum) Firmicutes strains further highlights a potential importance of metabolism of green-diet SQ by gut microbial communities to, ultimately, hydrogen sulfide.

Reconstructing the Evolutionary History of Chromosomal Races on Islands: A Genome-Wide Analysis of Natural House Mouse Populations.

Chromosomal evolution is widely considered to be an important driver of speciation, as karyotypic reorganization can bring about the establishment of reproductive barriers between incipient species. One textbook example for genetic mechanisms of speciation are large-scale chromosomal rearrangements such as Robertsonian (Rb) fusions, a common class of structural variants that can drastically change the recombination landscape by suppressing crossing-over and influence gene expression by altering regulatory networks. Here, we explore the population structure and demographic patterns of a well-known house mouse Rb system in the Aeolian archipelago in Southern Italy using genome-wide data. By analyzing chromosomal regions characterized by different levels of recombination, we trace the evolutionary history of a set of Rb chromosomes occurring in different geographical locations and test whether chromosomal fusions have a single shared origin or occurred multiple times. Using a combination of phylogenetic and population genetic approaches, we find support for multiple, independent origins of three focal Rb chromosomes. The elucidation of the demographic patterns of the mouse populations within the Aeolian archipelago shows that an interplay between fixation of newly formed Rb chromosomes and hybridization events has contributed to shaping their current karyotypic distribution. Overall, our results illustrate that chromosome structure is much more dynamic than anticipated and emphasize the importance of large-scale chromosomal translocations in speciation.

Diving into divergence: Differentiation in swimming performances, physiology and gene expression between locally-adapted sympatric cichlid fishes.

Sympatric speciation occurs without geographical barriers and is thought to often be driven by ecological specialization of individuals that eventually diverge genetically and phenotypically. Distinct morphologies between sympatric populations occupying different niches have been interpreted as such differentiating adaptive phenotypes, yet differences in performance and thus likely adaptiveness between them were rarely tested. Here, we investigated if divergent body shapes of two sympatric crater lake cichlid species from Nicaragua, one being a shore-associated (benthic) species while the other prefers the open water zones (limnetic), affect cruising (Ucrit ) and sprinting (Usprint ) swimming abilities - performances particularly relevant to their respective lifestyles. Furthermore, we investigated species differences in oxygen consumption (MO2 ) across different swimming speeds and compare gene expression in gills and white muscle at rest and during exercise. We found a superior cruising ability in the limnetic Amphilophus zaliosus compared to the benthic Amphilophus astorquii, while sprinting was not different, suggesting that their distinct morphologies affect swimming performance. Increased cruising swimming ability in A. zaliosus was linked to a higher oxygen demand during activity (but not rest), indicating different metabolic rates during exercise - a hypothesis supported by coinciding gene expression patterns of gill transcriptomes. We identified differentially expressed genes linked to swimming physiology, regulation of swimming behaviour and oxygen intake. A combination of physiological and morphological differences may thus underlie adaptations to these species' distinct niches. This complex ecological specialization probably resulted in morphological and physiological trade-offs that contributed to the rapid establishment and maintenance of divergence with gene flow.

Conservation and novelty in the microRNA genomic landscape of hyperdiverse cichlid fishes.

MicroRNAs (miRNAs) play crucial roles in the post-transcriptional control of messenger RNA (mRNA). These miRNA-mRNA regulatory networks are present in nearly all organisms and contribute to development, phenotypic divergence, and speciation. To examine the miRNA landscape of cichlid fishes, one of the most species-rich families of vertebrates, we profiled the expression of both miRNA and mRNA in a diverse set of cichlid lineages. Among these, we found that conserved miRNAs differ from recently arisen miRNAs (i.e. lineage specific) in average expression levels, number of target sites, sequence variability, and physical clustering patterns in the genome. Furthermore, conserved miRNA target sites tend to be enriched at the 5' end of protein-coding gene 3' UTRs. Consistent with the presumed regulatory role of miRNAs, we detected more negative correlations between the expression of miRNA-mRNA functional pairs than in random pairings. Finally, we provide evidence that novel miRNA targets sites are enriched in genes involved in protein synthesis pathways. Our results show how conserved and evolutionarily novel miRNAs differ in their contribution to the genomic landscape and highlight their particular evolutionary roles in the adaptive diversification of cichlids.